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Abstract 

We study the lifetimes of TeV-scale heavy neutral leptons (Majorana neutrinos) that appear 



in a model suggested by Okamura et al. 



a. 



We develop a convenient way to parametrize the 
neutrino mass texture of the model, and illustrate our method by calculating the mass spectrum, 
decay widths, and lifetimes of the heavy particles over the entire parameter space. From the mass 
spectrum, we find that for most of the parameter space, only two-body decays are relevant in the 
calculation of the lifetime, with typical values falling in the range of 10 -26 to 10~ seconds. If the 
particles discussed here are created at colliders, their lifetimes are short enough for them to decay 
inside the detector, while long enough to lead to a narrow peak in the invariant mass spectrum of 
the decay products. However, an analysis by Dicus, Karatas, and Roy suggests that they may 
be difficult to observe at the LHC. 
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I. INTRODUCTION. 



Several models of neutrino mass have been suggested in the literature in which the neu- 
trinos acquire masses through a seesaw type mass texture, but the Majorana masses of 
be right-handed neutrinos are at the TeV scale instead of the GUT scale of ~ 10 16 GeV 
. The smallness of the neutrino masses in those models is achieved either by the 
reduction of the rank of the mass matrix through a judicious choice of mass texture P, [J, 
or by the suppression of the Dirac masses through an extended Higgs sector 0,0,0]. 

In such models, the heavy, mostly-right-handed mass eigenstates typically have masses of 
a few TeV, placing them within reach of the CERN Large Hadron Collider (LHC) or future 
e + e~ linear colliders. If created, the particles will decay into a light neutrino+Higgs through 
the Yukawa interaction responsible for the Dirac masses, or into a light neutrino+Z or a 
charged lepton+VT through the small admixture of the left-handed neutrino state. This last 
decay mode is particularly interesting since the decay products can be all visible. Of course, 
whether such a decay, and thus the particle, can be observed at colliders or not depends on 
whether the lifetime of the particle is short enough for it to decay inside the detector, and 
if that is the case, whether the width is small enough so that a narrow peak is discernible 
in the invariant mass of its decay products. 

In this paper, we calculate the lifetimes of the heavy, mostly-right-handed states of the 
model proposed by Okamura et al. in Ref. Q]. The original motivation of the model was 
to explain the NuTeV anomaly HQ, one possible solution of which requires >arg 18 h mixing 
(6 2 ~ 0.003) between the light and heavy (3> M z ) neutrino states [9|. Denoting the left- 
and right-handed neutrino states by v and £, respectively, the Okamura texture is given by 
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where the dimensionless parameters a, (3, and 7 are in general complex and assumed to 
satisfy the relation 

a + p + j = 0. (2) 



This condition reduces the rank of the above mass matrix to three, leading automatically to 
three massless neutrino states. Though the actual light, mostly left-handed neutrino states 
in nature are not completely massless, this model suffices as a first approximation. We fix 
the normalization of the three complex parameters a, (3, and 7 to 



Ictl 2 + 



| 2 + l7| 2 



3 . 



(3) 



The dimensionful parameters m and M can be taken to be real and they set the scale of 
the Dirac and Majorana masses, respectively. The solution to the NuTeV anomaly requires 
their ratio to be P, 0| 



m 
M 



0.03 . 



(4) 



If the gauge singlet states & (i = 1, 2, 3) couple to other particles only through the Yukawa 
interactions which generate the Dirac submatrix of Eq. (fljl, then any permutation of the 
three complex parameters a, (3 and 7 leads to the exact same model since we will have the 
freedom to relabel the three gauge singlet states without affecting any physics. In those 
cases, there exist a 3! = 6 fold redundancy in the parameter space spanned by a, (3, and 7. 
This will be assumed in the following. 
If we set M = in Eq. we obtain 

am (3m 7m 
am (3m 7m 
am (3m 7m 



am am am 
(3m (3m (3m 
7m 7m 7m 

which is manifestly rank 2. The non-zero eigenvalues of this matrix are 1 



(5) 



imv/3 (|a| 2 + |/3| 2 + I7I 2 ) = ±3 m . 



(6) 



Therefore, this mass texture leads to four massless and two massive Majorana fermions. 
Pairing up the Majorana fermions with the same mass and opposite CP, we can reduce the 
set to one massive and two massless Dirac fermions jll3]. If we assume that the up-type 



1 A factor of \/3 is missing from Eq. (65) of Rcf. 2J. 
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quarks share the same Dirac mass texture as the neutrinos, as would be the case in the 



Pati-Salam model 11 1, we obtain one massive quark which can be identified with the t, and 
two massless quarks which can be identified with the u and the c. To produce the t quark 
mass, we need 

m ~ 60 GeV , (7) 

which together with Eq. (J1J) implies 

M ~ 2 TeV . (8) 

Fixing m and M to these values, the parameter space of the Okamura model is given by the 
values of a, f3, and 7 which satisfy Eqs. (0) and ©• 

In the following, we introduce a convenient graphical representation of the parameter 
space for the Okamura model, and calculate the masses and lifetimes of the three heavy 
mass eigenstates over it. We find that except for the vicinity of three isolated points at 
the 'edge' of the parameter space, the three masses are always in the TeV range, and the 
lifetimes are typically in the range of 10~ 26 to 10 -24 seconds. In terms of the widths, these 
correspond to the range of 0.7 ~ 70 GeV, which are fairly narrow compared to the masses. 



II. THE PARAMETER SPACE OF THE OKAMURA MODEL 

We begin by noting that for the three complex parameters a, /3, and 7 to sum to zero, 
Eq. (J2J), they must form a closed triangle when summed tip-to-tail as vectors in the complex 
plane. Without loss of generality, we can set the phase of a to zero. This can always be 
achieved by changing the overall phase of a, [3 and 7, and does not affect any physical result. 
Therefore, the triangle formed by a, /3, and 7 can be assumed to have its base along the 
positive real axis. We define the "orientation" of this triangle as the direction of the vectorial 
cross product a x (3. If the orientation of the triangle is (out of the complex plane), then 
(3 is in the upper complex plane while 7 is in the lower complex plane. If the orientation of 
the triangle is ® (into the complex plane), then (3 is in the lower complex plane while 7 is in 
the upper complex plane (see Fig. [TJ. Then, it is easy to see that specifying the lengths of 
the three sides \a\, \/3\, and I7I, and the orientation of the triangle is equivalent to specifying 
the three complex numbers a, (3, and 7. 



4 



Imz 



Imz 



Imz 












a 


— *• ► 


7 


t 






7 /\ 

/ 


R 

i 




a 






t 






/ 






\a\ 





Rez 




a ' 




% 




Rez 




equivalent 



a 



► Rez 





FIG. 1: The three complex numbers a, (5 and 7 satisfying a + (3 + 7 = form a close triangle. For 
each choice of the three lengths \a\, \/3\, and I7I, there are two possible orientations of the triangle 
(0 and 0) which are related by complex conjugation (1st and 2nd columns). However, the case 
is equivalent to the case with the lengths |/3| and |7| interchanged (2nd and 3rd columns). 



Furthermore, we need not consider both orientations since the two cases can be trans- 
formed into each other by a simple interchange of the lengths of (3 and 7, and a relabeling of 
the singlet neutrino fields. As discussed previously, this does not affect any physical result 
either. Therefore, we will always take the triangle to be in the orientation. This choice 
also reduces the redundancy of the parameter space from 3! = 6 to 3 since we have used up 
the freedom to interchange j3 and 7 to fix the orientation. 

This consideration shows that specifying the three lengths |a|, \/3\, and I7I suffices to 
uniquely determine the Okamura texture, with cyclic permutations of the three lengths 
leading to the same model. (This residual redundancy comes from our freedom to choose 
which of the three lengths to call \a\.) The question then, is, how can we specify those three 
lengths so they satisfy the normalization condition Eq. ©, and also the triangle inequalities: 
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FIG. 2: The distances from any point inside an equilateral triangle of height three to its three sides 
add up to three. We can use these distances to specify \a\ 2 , \(3\ 2 , and I7I 2 . The triangle inequality 
is satisfied for points inside the unit circle which inscribes the triangle. 



\P\ + \j\>\a\, \a\ + |/3| > | 7 | , |a| + | 7 | > \P\ , (9) 

so they form a closed triangle? To this end, we utilize the fact that the sum of distances 
from any point inside a triangle to its three sides is constant: any point inside an equilateral 
triangle of height three will have distances to the three sides which add up to three. If 
we identify these distances with \a\ 2 , \/3\ 2 , and I7I 2 , we can use the position of the point 
to specify the three lengths. Requiring the square-roots of these distances to satisfy the 
triangle inequality constrains the point to be inside a unit circle which inscribes the triangle. 
Therefore, for every point inside the unit circle, we can associate a corresponding parameter 
set for the Okamura texture (see Fig. |2J). 

If we specify the position of a point inside the unit circle with its polar coordinate (r, 9) , 
where < r < 1, 6 G [—71,71), the corresponding values of \a\, |/3|, and | 7 | are: 



a 




Vl + r sin 9 , \/3\ = Jl + rsin ( 6- — ) , | 7 | = x jl + r sin ( 6 + — J . (10) 
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The phases of the three numbers are: 



arg a 
arg (3 

arg7 







7T — COS 



1/2 + rsin(fl - vr/3) 



-7T + COS 



y/[l + rsin#][l + rsin(9 - 2vr/3) 
l/2 + rsin(# + 7r/3) 



, = . (11) 

y/[ 1 + r sin 9 ] [ 1 + r sin(# + 2tt/3) ] 

In Fig. El we plot the dependence of arg (3 and arg 7 on the position of the point inside the 
unit circle. 

A cyclic permutation of a, (3, and 7 which leaves the physics invariant up to an overall 
phase corresponds to the transformation 9 — > 9 + 27r/3 (120° rotations). This means that 
we expect the same symmetry to be present in the mass spectrum and the values of heavy 
neutrino decay widths and lifetimes. This can be used as a useful check of our calculations. 




■ 




FIG. 3: Contour lines and density plots for arg/3 (left) and arg 7 (middle). The distances between 
two consecutive equipotential lines are Aarg/3 = A arg 7 = 0.057T. The color scheme is shown on 
the right. 



III. THE LAGRANGIAN 



To calculate the lifetimes of the heavy neutral states, we must first specify their interac- 
tions. We denote the left-handed charged lepton fields with £, and the left- and right-handed 
neutrino fields with v and £, respectively: 



£2 

6 



(12) 
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The right-handed neutrino fields, & (i — 1, 2, 3), are gauge singlets. The components of the 
Higgs doublet are denoted 



H 



00 



(13) 



Then, the Lagrangian which governs the interaction of the neutrinos is 



C — Cw,z + £h + £m , 



(14) 



where 
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9 m 
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c H = -z x ((f) v - <p + e) + h.c. , 

C M = -\lMi c + h.c. . 



2 cos 9 W 



(15) 



We neglect the Yukawa interactions which give rise to the charged lepton masses: the charged 
leptons are treated as massless as well as the light neutrino states. In the Okamura model, 
the Yukawa matrix A and the Majorana mass matrix M. are given by 



A = 



V2- 



rn 
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M = M 
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(16) 



After the neutral Higgs develops a VEV, 



y/2' 



(17) 



the Yukawa matrix A leads to the Dirac mass matrix of the neutrinos: 



V = —=X = m 
V2 
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(18) 



The Goldstone bosons are absorbed into the W and the Z, 0° — > 1 / \/2(h + v) as usual, and 
the resulting Lagrangian is: 
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-ZVu--j=(t\u)h-±tM? + h.c 



(19) 



The neutrino mass terms can be written as 



ZVv + -ZMZ c + h.c. 



-{^Vu + ^V T C + iMC)+h.c. 
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(20) 



This mass matrix is diagonalized with a unitary transformation involving the v and £ c fields: 

(21) 
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(22) 



with M diag = diag(0, 0, 0, M 1 , M 2 , M 3 ). The rj and x fields are the left-handed mass eigen- 
fields with rj being the light (massless) states, and x being the heavy states. Decomposing 
the 6x6 matrix U into four 3x3 matrices as 



U 



(23) 



we can write 



v = U uv r] + U vx x , 

e = u^ v c + ui xX c , 



(24) 



(Because the r\ fields are exactly massless and degenerate in our model, the matrices U m) 
and are not uniquely determined. However, this does not affect our final results. Note 
also, that though U is unitary, its four 3x3 submatrices are non-unitary in general.) The 
relevant interaction terms in the Lagrangian involving the x fields are then: 
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2 cos 9 W 
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viut^Yx + xiu^u^Yv z„, 
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9 m 
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— (£\u)h^ -^[7f{ulXU ux )x + r(UlXU uv ) V \h, (25) 
plus the Hermitian conjugates of the later two lines. Introducing the Majorana fields 



n — rj + rf 



N = x + X c 



(note that these fields do not have definite lepton number) we can write 

V = P L n, v c = P R n, X = P L N } X ° = PrN , 

and the relevant interaction Lagrangian in terms of these fields becomes 

C = —^-\n{AfP L - A*YPr)n}z, 
2 cos W i J 

+ (1B^P L N) W~ - -j= (FB*YPrN) W; 
-n(chP L + C*~hP R ) N , 



where 



B = U, 



We have used the generic relations [3| 

^1 OP RtL ^2=^2 T P R , L ft , ?1 O^P R)L ^2 
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(27) 



(28) 



(29) 



(30) 



(O is a matrix which carries flavor indices only), and the fact that n c = n and N c = N by 
construction, to rearrange the terms in Eq. (}2*Hj) in such a way that all the iV-fields stand 
at the rightmost position of each term to facilitate the extraction of the iV-decay matrix 
elements. 



IV. LIFETIMES 

From Eq. (}2*8*j) . we can immediately derive the amplitudes for the two-body decay pro- 
cessed of the heavy neutrinos, iVj (i = 1, 2, 3), shown in Fig. 0] If the iVj were lighter than 
the W, Z, or h, then we will need to consider three-body decay processes mediated by these 
particles, but it turns out that except for small neighborhoods around isolated points in the 
parameter space of the model, they are always heavier. It therefore suffices to consider only 
the two-body decay modes. 
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FIG. 4: The 2-body decay processes of the heavy neutrino N. 

Now, straightforward calculations allow us to write down the partial decay widths for 
each channel of decay (the indices % and j below run from 1 to 3): 
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16tt " V" M i 

The first two lines can be compared with the results of Djouadi in Ref. [13J. At first sight, 
these expressions may seem to imply that the iV — > nh channel is suppressed with respect to 
the other two since its partial width grows linearly with the mass Mj, while for the N — > nZ 
and iV — > IW channels the widths grow as Mf. However, the interactions of the heavy 
Majorana neutrino N with the gauge bosons are suppressed because only a small fraction 
of N is the left-handed neutrino v. Since most of N is the right-handed neutrino £, no such 
suppression exists for its interaction with the Higgs h. Numerically, it turns out that all 
three channels of decay must be taken into account. 



V. RESULTS 



Now we have everything at hand to calculate the lifetimes of the heavy neutrinos in the 
Okamura model. The parameter space of the model is represented by the interior of a unit 
circle as discussed in section |H] For each point inside the unit circle, we can calculate the 
Okamura texture using Eqs. ((TJ), (J7J), (jSJ), (jTUjl . and ([lip, diagonalize it to obtain the masses 
and mixings Q| , and calculate the decay widths and lifetimes of the heavy neutrinos using 
Eq. pip. As the Higgs mass, we take m# = 200 GeV. (The choice of the Higgs mass has 
little effect on our result as long as m H < Mj.) 
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FIG. 5: (a), (b), (c) density and contour plots for masses Mi, M2, and M3 of the lightest N\, 
medium heavy N2 and heaviest N$ heavy neutrino respectively (TeV). The distances between two 
consecutive equipotential lines are AMi = 0.1 TeV, AM2 = AM3 = 0.05 TeV; (d) the vicinity of 
the point where Mi approaches zero (r = 1, 9 = —tt/2); A Mi = 0.025 TeV; (e) mass color coding. 

The resulting contour and density plots for masses, decay widths, and lifetimes of the 
heavy neutrinos iVj, N 2 and N 3 are presented in Figs. EH3 First, note that the graphs are 
symmetric under rotations by multiples of 27r/3 as anticipated in section |H] Next, from 
Fig. El we can easily see that the values of the heavy neutrino masses are larger than the W, 
Z, or Higgs thresholds for most of the parameter space, justifying our use of two-body decay 
amplitudes. The mass of N\ becomes smaller than these thresholds only in the vicinity of 
three isolated points at r = 1, 6 — —n/2 + 2nk/3 (k = 0, 1, 2), as illustrated in Fig. 01. As 
was shown in Ref. Q], at these three points one of the iV-nelds is completely massless while 
the other two have degenerate mass. The lightest N particle is completely stable at these 
points with zero decay width and infinite lifetime. However, the existence of such a lig ht 
(less than W and Z thresholds) N particle is already ruled out experimentally by L3 jl5( so 
we need not consider these points further. 
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FIG. 6: (a), (b), (c) density and contour plots for widths Ti, T2, and T3 of the lightest JVi, 
medium heavy N2 and heaviest N3 heavy neutrino respectively (GeV). The distances between two 
consecutive equipotential lines are AFi = 1 GeV, = 2 GeV, = 3 GeV; (d) the detailed 
picture of the central part of T\; AT\ = 1 GeV; (e) width color coding. 

It was also shown in Ref. [jj] that at the center of the circle one heavy neutrino completely 
decouples from the light neutrino states (and therefore from the rest of the Standard Model 
particles) while the other two heavy states have degenerate masses. This decoupling can be 
seen in Fig. Eli and Fig. [7|i where at the center of the circle the decay width of the lightest 
heavy neutrino is zero and the lifetime is infinite. A similar decoupling occurs at the points 
where r = 1 and 9 = —5tt/6 + 27r/c/3, k — 0, 1, 2. 

Except for the vicinity of these points, the lifetimes of the N particles are typically in the 
range of 10 -26 to 10 -24 seconds (see Fig. |7|). Assuming that the particles are non-relativistic, 
the maximum distance they can travel from their production points before decay is in the 
range of 1CT 17 to 10~ 15 meters. If produced at colliders, they will decay inside the detector. 
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FIG. 7: (a), (b), (c) density and contour plots for lifetimes tl, T2, and T3 of the lightest N±, 
medium heavy N2 and heaviest N3 heavy neutrino respectively (10 -26 sec). The distances between 
two consecutive equipotential lines are An = 1 x 10~ 26 sec for t\ < 10 x 1CP 26 sec and for t\ > 
10 x 10~ 26 sec equipotential lines corresponding to t\ = 20 x 10~ 26 , 40 x 10~ 26 and 200 x 10~ 26 sec 
are presented, At2 = 0.2 x 10~ 26 sec, Ar3 = 0.1 x 10~ 26 sec; (d) the detailed picture of the 
central part of n> Arx = 10 x 10~ sec for n < 100 x 10" 26 sec and An = 100 x 10" 26 sec for 
n > 100 x 10~ 26 sec; (e) lifetime color coding. 

On the other hand, the width-to-mass ratios of the particles are in the range of 0.1 to 3 
percent as shown in Fig. |HJ Therefore, the invariant mass spectrum of the decay products 
can be expected to show a very narrow peak. 
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FIG. 8: (a), (b), (c) density and contour plots for the widths-to-mass ratios of the lightest N%, 
medium heavy N2 and heaviest N3 heavy neutrinos respectively. The distances between two con- 
secutive equipotential lines is 0.001 for all plots; (d) mass-to-width ratio color coding. 

VI. SUMMARY AND DISCUSSION 



In this paper, we calculate the mass spectrum, decay widths and lifetimes of the mostly- 
right-handed heavy neutrino states that appear in the model proposed by Okamura et al. 
in Ref. Q|. We map the parameter space of the Okamura model to the interior of a unit 
circle, and represent the results of our calculations as density-contour plots over it. For 
the phenomenologically viable region of the model's parameter space, the heavy states have 
masses of a few TeV, and are short-lived with the typical lifetimes from 10~ 26 to 10~ 24 
seconds. At the same time, the decay widths are very small comparing to the masses. The 
typical width-to-mass ratio is in the range of 0.1 to 3 percent. 

Though we have found that the heavy neutral particles in the Okamura model have 
lifetimes in the r ange that allow for their observation at colliders, an analysis by Dicus, 



Karat as, and Roy 



l(j| suggests that they may be difficult to observe at the LHC. In Ref. 



[16], 



the authors consider the production of like-sign leptons, a lepton number violating process, 
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as the signature of the heavy mostly-right-handed Majorana neutrino N: this can occur 
through the process — > £^N — > £ ± £ ± W /=F — > £ ± £ ± + jets, or the t- and w-channel 
exchange of A" between two like-sign W's radiated from the protons. Assuming masses in 
the range of 150 ~ 4000 GeV, and a mixing as large as 6 2 ~ 0.043 (an order of magnitude 
larger than the Okamura model), the number of expected events at the LHC was shown 
to be only a few per year; too few for them to be discernible above the Standard Model 
background. 

If the gauge group is extended to SU{2) L x SU{2) R x U(1)b-l, then the A^'s can be 
copiously produced through the Wr and Z', as discussed in Refs. 0, Q, Q- We will 
discuss the embedding of the Okamura texture into this gauge structure in a subsequent 
paper j^j]. 
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